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Abstract 

Damping of an electromagnetic wave in a strong magnetic field is an- 
alyzed in the kinematic region near the threshold of electron-positron pair 



Q I production. Damping of the electromagnetic field is shown to be noticeably 

1-^ ' nonexponential in this region. The resulting width of the photon 7 — > e~^e~ 

decay is considerably smaller than previously known results. 
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The problem of propagation of electromagnetic fields through an active medium 
is inherent in a variety of physical phenomena. The birth and evolution of supernova 
and neutron stars, where the matter density can be on the order of nuclear density 
p ~ 10^^ — IQ^^g/crn? and the temperature can achieve several tens of MeVs, are 
the largest scale and the most interesting such phenomena. In addition to dense 
and hot matter, a strong magnetic field, which can be several orders of magnitude 
as high as so-called critical, or Schwinger, value B^. = m^/e ~ 4.41 ■ lO^^G Q, can 
be generated in the above-mentioned objects [^ 0]. This strong magnetic field 
can induce new phenomena which can considerably affect the evolution of these 
astrophysical objects. Electromagnetic-field damping caused by electron-positron 
pair production in an external magnetic field is one of these phenomena. Recall 
that the 7 -^ e^e~ process is kinematically forbidden in vacuum. The magnetic 
field changes the kinematics of charged particles, electrons and positrons, allowing 
the production of an electron-positron pair in the kinematic region qf^ = q^ ~ ^3 — 
4ml, where go is the photon energy and ^3 is the momentum component along the 
magnetic field | 

In 1954, Klepikov [Q examined the production of an electron-positron pair by a 
photon in a magnetic field and obtained the amplitude and width of the 7 —>■ e^e~ 
decay in the semiclassical approximation. Later, the authors of 0-0 considered 
this process in the context of its astrophysical applications. It was pointed out 
in [0, ^ that the use of the expression derived in for the width considerably 
overestimates the result in the strong magnetic field limit. In this case, one should 
use an exact expression for the width of one-photon production of a pair when 
electrons and positrons occupy only the ground Landau level. However, it was 
found that the expression for the decay width in the limit of strong magnetic field 
has a root singularity at the point gj^ = 4 ml. Shabad [^] emphasized that 
this behavior indicates that the decay width calculated in the perturbation theory 
cannot be treated as a damping coefficient. In this case, the damping coefficient is 
primarily determined from the time evolution of the photon wave function in the 
presence of a magnetic field. Shabad suggested that this dependence be obtained 
by solving the dispersion equation with account taken of the vacuum polarization 
in a magnetic field with complex values of photon energy. In our opinion, this 
method has several disadvantages. First, it is well known but rarely mentioned that 
the dispersion equations with complex energies have no solutions in the physical 
sheet. Solutions are in the nonphysical Riemannian sheets (analyticity region of 

^We use the system of units where Ti = c— 1. 

^Hereafter, we consider the magnetic field directed along the third axis. 



the polarization operator), which are generally infinite in number. As a result, 
the dispersion equation has the infinite number of solutions with both positive and 
negative imaginary parts of energy. The physical status of these solutions requires 
a separate investigation. 

Shabad |§ used the asymptotic form of the polarization operator near the pair 
production threshold and erroneously treated it as a two-sheet complex function. 
This circumstance led to the existence of two complex conjugate solutions, one of 
which is physically meaningless because it has a positive imaginary part and, there- 
fore, provides exponentially increasing amplitude of electromagnetic wave. There- 
fore, to obtain physically meaningful result, one should artificially discard the re- 
dundant solutions. 

Second, this approach cannot correctly describe the substantially nonexponential 
damping in the near-threshold region in a strong field. 

Thus, damped electromagnetic waves in a magnetic field cannot be completely 
described by solving the dispersion equation. 

In this work, we use a method that is applied in the field theory at finite tem- 
peratures and in plasma physics (see, e.g., [|10|)- It consists of the determination of 
a retarded solution to the electromagnetic field equation that includes an external 
source and takes into account the vacuum polarization in a magnetic field. Time 
damping of the electromagnetic wave is analyzed in a uniform external magnetic 
field, whose intensity is the largest parameter of the problem, B^ ^ q'^,ml. To 
describe the time evolution of electromagnetic wave Aa{x) [x^ = (t,x)], we consider 
a linear response of the system {Aa{x) and a vacuum polarized in magnetic field) 
to an external source, which is adiabatically turned on at t = — cxo and turned off at 
t = 0. At t > 0, the electromagnetic wave evolves independently. Thus, the source 
is necessary for creating an initial state. For simplicity, we consider the evolution 
of a monochromatic wave. In this case, the source function should be taken in the 
form 

Ja{x)=jae''^e''9i-t), £^0+. (1) 

The time dependence of Aa{x) is determined by the equation 

{gap dl - d^dp) Ap{x) + j d^x' Va^ix - x') Ap{x') = J^{x), (2) 

where Vap{x — x') is the photon polarization operator in a magnetic field. We note 
that, for the source on the right-hand side of Eq. (Q) to be conserved, daJa = 0, 
the current ja must have the form ja = (0, j), j ■ k = 0. The evolution of Aa{x) is 



described by the retarded solution 

A^{x)= l<fx'G%{x~x')Mx'), (3) 



where G^Jx — x') is the retarded Green's function, which is defined through the 
commutator of the Heisenberg operators of electromagnetic field as (see, e.g., |11|) 

Gf^(x-x') = -z(0|K(x), Ap{x')]\0)9it-t'), (4) 

It is instructive to express this function in terms of the causal Green's function 

C^pix-x') = -z(0|Ti,(x)i^(x')|0), (5) 

by using the relationship 

Gf^(x - x') = 2ReG'^f,ix - x') e{t - t'). (6) 

In the presence of a magnetic field, it is convenient to decompose Green's function 
(Bl) in the eigenvectors of polarization operator [H : 



^^M^) = jAiGU'i)^-'''' (7) 
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where V'^^^q) are the eigenvalues of polarization operator. The eigenvectors 

b^a'> = (q^U (9) 

ba^ = q^{qVV)a - {q^^q)qa, (10) 

together with the 4-vector g„ form a complete orthogonal basis in the Minkowski 
4-space. In Eqs. ([T0|) , yj^/? = Fajs/B, (p^p = ^Sais^iuf^iu are dimensionless magnetic- 
field tensor and dual tensor, respectively, {q^)a = qa^^aa, {qwq) = qa^ap^iiaqa- 
Substituting Eqs. (|l]) and (H) into Eq. (|) and using Eqs. (|^) and (H), we obtain 
after simple integration the following result at t > 0: 

4 



where qa = (i?o;k). Note that the definition of the integral in Eq. ( [TT| ) should 
be completed because the integrand can include singularities, which are due, on 
the one hand, to zeros of its denominator and, on the other, to the domain of its 
definition. To analyze these singularities, it is necessary to know the explicit form 
and analytical properties of the eigenvalues V^'^\q) of the polarization operator, 
which was examined in detail in a number of works. In the limit of strong magnetic 
field, the functions V^'^^q), which we are interested in, can be borrowed, e.g., from 



T^, |T^ and represented as [with the 0(1/ B) accuracy] 



P«(g) 



a 



-^ql-q^A(B,q'), 
on 

q'AiB,q'), 



2a eB 



TT 



z(l-z) 



: arctan , 



q'A{B,q'), 



(12) 
(13) 
(14) 



where 



A{B,q' 



[1.792 - \n{B/B,)] + n{q'^ 



2 I 2 
?1 +?2) 



a 
3^ 
z = ql/Aml, ql 

and q^'Tc{q^) is the photon polarization operator in the absence of a magnetic field. 
Note that the contribution from the pole (jj^ = that results from the normalization 
of the basis vectors h^^^ and fo^^-* is nonphysical and, taking into account explicit form 
(P^)-(|T^) of the polarization operator, can be removed by gauge transformation after 
summation over polarizations. Thus, the contribution to the solution can be made 
only by the poles corresponding to the dispersion equation 



g2-p(^)(g) = 0. 



(15) 



Using solution 
two modes, A = 



l^pl]), one can demonstrate on the basis of Eqs. (|T^)- ([l^) that only 
1 and A = 2 with the polarization vectors 



r(l) 



5(1) 



{(iV)a, 



r(2) 



&(^) 



(#)a 



[hW] 



(U^yy 



are physically meaningful for real photons f\ 



■^ Modes with the polarizations Sa and Eq correspond to so-called parallel 



(16) 



I ) and perpendic- 



ular (_L) modes, respectively, in the Adler notation ||14 



A photon of the third, A = 3, mode is nonphysical 0. Indeed, substitution 
of the expression for V^^\q) into Eq. (^) gives the equation that has the only 
solution q^ = 0. Therefore, the contribution of the third mode to solution (^) is 
proportional to the total divergence and can be eliminated by the corresponding 
gauge transformation. 

In the limit of strong magnetic field, only the mode with the polarization vector 
e^^^ can decay into an electron-positron pair, because only the eigenvalue of the 
polarization operator V^'^\q) (|l^ has the imaginary part at g? > Ami- Therefore, 
to analyze time damping of the electromagnetic field, it is sufficient to consider only 
the term with V^^ (x) in Eq. ([TT|) . 

The further calculations can be simplified by going over to the reference frame, 
where k = {ki,k2,0), which can always be done without disturbing the properties 
of the external magnetic field. In this frame, g? = q^ and the polarization vector of 
the second mode takes the form e^^^ = (0, 0, 0, —1). As a result, V^'^\x) is expressed 
as 



where 



F{t) = I ^ 



-iqot 



C 



27rz {q,-^e){qi-k^-V(^\q)y 



V;(2)(0,x) = 242)^3e^kxj^g^^Q)_ 



The path of integration C in Eq. ([Tq) is determined by the analytical properties 
of V^'^\q) and is shown in Fig. 1. The function V^'^\q) is analytical in the complex 
plane go with cuts along the real axis (see Fig. 1). In the kinematic region |go| < 2me, 
the eigenvalue P*^^) (g) is real and Eq. ( jTSD has real solutions which govern the photon 
dispersion in this region. 

For further analysis, it is convenient to transform the path of integration to the 
path shown in Fig. 2. In this case, the integral in Eq. (llSD is represented as 



F{t)=Fpole{t)+F,ut{t), (19) 

where the first term is determined by the residue at the point go = cu, which is the 
solution to dispersion Eq. ([T5|). This term corresponds to the undamped solution 
in the region uj < 2?7ie H- The second term determines the time dependence of the 



electromagnetic field above the threshold of electron-positron pair production and 
has the form of the Fourier integral 

Fcutit) = j '^ F,^Mo)e-'''\ (20) 

— oo 

where 



R ^ ReV^'\qo) = ^eB( ^ ^ ln^+^^ + 2), (22) 



-1) v^-V^ 



Expressions (pOD -(|23D together with Eq. (|T^) determine the time evolution of the 
photon wave function above the pair production threshold in a strong magnetic field. 
Strictly speaking, because of the threshold behavior of the Fourier transform 
Fcutiqo), time damping of the function Fcutif) and, therefore, of the wave function 
Aa{t) is nonexponential. However, in some characteristic time interval (the inverse 
effective width of the 7 — >■ e^e~ decay can naturally be chosen as such an interval), 
the time dependence of the wave function can approximately be represented as 
exponentially damping harmonic oscillations: 

A^{t) ~ e-^=ff*/2 cos(cu.fft + 0o). (24) 

Here, oj^s and '-y^^ are, respectively, the effective frequency and width of photon de- 
cay, which should be found by using Eqs. (|2^)-(^B|) for each value of momentum 
k to determine the effective photon dispersion law above the threshold of electron- 
positron pair production. The quantity 7cff, which governs the intensity of photon 
absorption due to e"'"e~ pair production in a magnetic field, is important for astro- 
physical applications. The absorption coefficient obtained from the 7 — » e+e~ decay 
probability and containing a root singularity is usually employed in astrophysics 
(see, e.g., [0). Shabad [^ pointed out that this leads to the overestimation of 
the intensity of e"'"e~-pair production. Our analysis demonstrates that the calcu- 
lation of the absorption coefficient (decay width) by using the complex solution in 
the second Riemannian sheet also leads to a considerable overestimation in the 
near-threshold region, as is seen from Figs. 3 and 4. 



Nonexponential damping in the near-threshold region is known for the processes 
in vacuum and matter [^, |l^. However, as far as we know, it has not been consid- 
ered in an external field so far. In contrast to vacuum or medium, the near-threshold 
effect in the magnetic field is kinematically enhanced due to the singular behavior 
of the polarization operator in this field. Therefore, this phenomenon is not only 
topical for astrophysical application but is of conceptual interest. 
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Figure 1: The path of integration C in the complex go plane. The crosses are the 
poles corresponding to the real solutions of dispersion Eq. (|T3p. The shaded parts 
of the real axis are cuts. 
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Figure 2: The path of integration C after the transformation allowing one to separate 
the pole Fpoie{t) and cut Fcut{t) contributions. 
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Figure 3: The frequency dependence of the 7 -^ e~^e~ decay width in the near- 
threshold region for the magnetic field B = 200i?e- Line 1 is the tree approximation 
including the root singularity; line 2 is obtained from the complex solution of the 
dispersion equation in the second Riemannian sheet P[; and line 2 is 7eflf from 



approximation (^) 



0.8 



0.6 



0.4 



0.2 



100 200 300 400 500 

B/Be 

Figure 4: The decay width vs. the magnetic field for the frequency uj = 2.5me. The 
meaning of fines is tfie same as in Fig. 3. 



